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Abstract 

In the hot and dense hadronic sources formed in high-energy heavy-ion collisions, the particle in¬ 
teractions in medium might lead to a squeezed back-to-back correlation (BBC) of boson-antiboson 
pairs. We calculate the BBC functions of for sources evolving hydrodynamically in (2 -|- 1) 
dimensions and with longitudinal boost invariance. The BBC functions for hydrodynamic sources 
exhibit oscillations as a function of the particle momentum because the temporal distributions of 
hydrodynamic sources have sharp falls to 0 at large evolving times. The dependences of the BBC 
functions on the directions of the particle momentum are investigated. For transverse anisotropic 
sources, the BBC functions are minimum when the azimuthal angles of the particles reach 0. The 
BBC functions increase with decreasing absolute value of the particle pseudorapidity. The oscilla¬ 
tions and the dependences on the particle azimuthal angle and pseudorapidity are the significant 
signatures for detecting the BBC in high-energy heavy-ion collisions. 
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I. INTRODUCTION 


In the hot and dense hadronic sonrces formed in high energy heavy ion collisions, the 


mass modihcation of particles in mec 


(BBC) of the boson-antiboson pair 


0 


inm can lead to a sqneezed back-to-back correlation 


2|. This BBC is the result of a quantum mechanical 
transformation relating in-medium quasiparticles to the two-mode squeezed states of their 
free observable counterparts, through a Bogoliubov transformation between the creation 
(annihilation) operators of the quasiparticles and the free observable particles ll-l3|. The 
investigations of the BBC of boson-antiboson pairs may provide a new way for people to 
understand the thermal and dynamical properties of the hadronic sources in high energy 
heavy ion collisions. 

Denote Uk (aj^.) the annihilation (creation) operator of the free boson with momentum k 
and mass m, and 6k (6k) the annihilation (creation) operator of the corresponding quasipar¬ 
ticle with momentum k and modihed mass ttz* in homogeneous medium; they are related by 
the Bogoliubov transformation 0.01 


Ok 


= Ck 6k + s\ 6lj^, 


where 


Ck = cosh /k , Sk = sinh /k , /k = 2 log('^k/^k), 


cck = vk^ + nTg 


The BBC function is dehned as [l, 2| 

^(k, -k) = 1 + 




|a(k, -k)p 




a(k,k)G,(-k, -k)’ 

where Gc(ki,k 2 ) and Gs(ki,k 2 ) are the chaotic and squeezed amplitudes, respectively. 


( 1 ) 

( 2 ) 

(3) 

(4) 


Gc(ki,k2) = V^kitCk2 (aici^ka), 

G's(ki,k 2 ) = v'Dkphl^ (OkiOka), 


(5) 


( 6 ) 


where (■ ■ ■) indicates the ensemble average. The BBC 
with volume V and temperature T can be expressed as 


'unction for a homogeneous source 


C(k, -k) = 1 + 


c |ck 4 ^k + C-k s\ (n_k + 1) P 
V [ni(k)ni(-k)] 


(7) 
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where 


1 


rik = 


exp((]k/7') - 1’ 


ni(k) = |ckpnk + |s_kP(n_k + !)• 


( 8 ) 

(9) 


In Ref. j^, S. Padula et al. put forward the formulism of the BBC function for the local- 
equilibrium evolving system and studied the BBC functions of 00 for expanding sources 
with a Gaussian space prohle. Recently, the BBC functions of K^K~ were investigated J] 


for expanding Gaussian sources, and a method was suggested p to search for the squeezed 
BBG in heavy-ion collisions at the Relativistic Heavy Ion Gollider (RHIG) and the Large 
Hadron Gollider (LHG). In Refs. 6] and j^, we calculated the BBG functions of relativistic 
00 and K^K~ pairs for spherical and ellipsoid expanding Gaussian sources. The relativistic 
effect on the BBG functions 6(] and the dependence of the BBG functions on the direction 
of the particle momentum for the anisotropic sources ^ are investigated. However, all 
the source density distributions used in the calculations in previous works are space-time 
separated, having Gaussian spatial distributions and an independent temporal distribution 
of exponential decay or others js-^. Investigations of the BBG based on more realistic 
space-time evolving source models will be of interest in high-energy heavy-ion collisions. 

Relativistic hydrodynamics has been extensively applied to high-energy heavy-ion col¬ 
lisions. In this work, we use the ideal relativistic hydrodynamics in 2 -|- 1 dimensions to 
describe the transverse expansion of sources with zero net baryon density and combine the 
Bjorken boost-invariant hypothesis js] for the source longitudinal evolution. As a hrst step 
to the study of the BBG for a more realistic source model, these descriptions are suitable for 
the heavy-ion collisions at the RHIG top energy and the LHG energy j9l-l2l|. We investigate 
the BBG functions of 0 meson pairs for the hydrodynamic sources. The results indicate 
that the BBG functions of 00 exhibit oscillations as a function of the particle momentum 
and vary with the particle azimuthal angle and pseudorapidity. The oscillations and the 
dependences on particle azimuthal angle and pseudorapidity are the signihcant signatures 
for detecting the BBG in high-energy heavy-ion collisions. 

The rest of this paper is organized as follows. In Sec. H, we present the calculation 
formulas of the BBG function for hydrodynamic sources. In Sec. HI, we investigate the 
BBG functions of 00 for hydrodynamic sources with different initial geometries and energy 
densities. The oscillations of BBG functions as a function of the particle momentum and the 
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dependences of BBC fnnctions on the directions of the particle momentnm are also discussed 
in this section. Finally, a summary and conclusions of this paper are given in Sec. IV. 


II. CALCULATIONS OF BBC FUNCTIONS FOR HYDRODYNAMIC SOURCES 


The description of ideal hydrodynamics for the system with zero net-baryon density is 


dehned by the local conservations of energy and momentum 

d^T^^r) = 0 , 


10 |, 


( 10 ) 


where = [e{r)+V{r)] u^{r)u'^{r) —V{r)g'^^ is the density tensor of energy-momentum 

ideal fluid at space-time coordinate r, e(r) and V{r) are the energy density and pressure in 
the local rest frame of the fluid element at r, which move at velocity v{r), = 7 ( 1 , v) is 

the four-velocity, 7 = (1—and g^'' = diag(-|-, —, —, —) is the Minkowski metric tensor. 
Under the assumption of Bjorken longitudinal boost invariance js], we need only to solve 
the transverse equations of motion in the z = 0 plane, and the hydrodynamic solutions at 
z ^ 0 {v^ = z/t) can be obtained by the longitudinal boost invariance hypothesis QQ- 
From Eq. (ITOl) we have the transverse equations in the z = 0 plane. 


( 11 ) 


dt£ + d,[iS+V)v^]+dy[iS+V)vy] = 
dtM^ + d,iM^v^+V) + dy(M^v^) = -Q{M\ t), 
dtM> + V) + dy{M>v^\V) = -Q{M\ t), 

where ^ = 'j‘^{e + V)-V, = -/^{e + Vy, {i = x,y), t) = {S + V)/t, 

and Q{A4\t) = M.yt. In equation set (ITTll there are e, P, and four variables. So 
an equation of state, P(e), is needed to enclose the equation set. In the calculations, we 
use the equation of state of s95p-PCE, which combines the hadron resonance gas at low 


temperatures and the lattice QCD results at high temperatures 


20l |. We assume that the 


system reaches the static local equilibrium at tq = 0.6 fm/c after the collision, and take the 
initial energy density distribution in the transverse plane as the Gaussian distribution. 


e = e„exp[-x^/(2R^) - !/ 7 ( 2 fiJ)|, 


( 12 ) 


where Cq and Ri (i = x, y) are the parameters of the initial source energy density and radii. 
With the equation of state and the initial energy density we can solve equation set .CD 


using the relativistic HLLE scheme and Sod’s operation splitting method [9, 


13|, l2ll-l26| 
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For hydrodynamic sources, with the formula derived 
chaotic and squeezed amplitudes can be expressed as 

G'e(ki,k 2 ) = 


w Makhlin and Sinyukov [27[, the 

3| 



* ^iqi, 2 -rj 

\J 1 

1 (27r)3 

^1,2^ ] 

Tkpk'l 


^n' 


V v 

Ki,K2 


I S 1 ^/ u/ I [ U 1 ^/ 1 ^/ 

I K.2 I L K.2 


+ 


(13) 


G.(ki,k 2 ) = 


(27r)3 


-'^ 1,2 


2i K\ 2-fJ 


X77- ^\c' \c' ^ \c' 

Ki,K2 k.^, K2 


n 


V' —l<r' 


+ 


u}. 


( 14 ) 


Here d^a^{r) is the four-dimension element of freeze-out hypersurface, qi 2 = k^—ki^, K ’^2 = 
(fcf -|- ^ 2 ) 72 , and k' is the local-frame momentum corresponding to kj {i = 1, 2). The other 
local variables are: 


c±k',±k' = cosh[/ 4 ki ], 


/±kp±k'2 “ 2 
1 


■^ikpik^ — sinh[/4^^±k^ ], 
1 


(wh +ujL)/{Q!]^, +^k') 


= 2 [Kl2U^.{r)/Kl'^2Uu{r)\ 

= /ki,k 2 (r). 


(15) 

(16) 


(17) 


oo^{r) = -^kf (r) + m? = k'^'u^^r) 
= 7 v[^^ki -k, ■ v(r)], 

n'k'(r) = 

= \J[kiU^{r)Y — w? + 1711 
= kru,{r), 


(18) 


(19) 


n. 


±kp±k' = exp < - 


2 (^k; + “ hl, 2 (^) 

= exp {-[iFi7>;,(r) - /ii, 2 (r)] / T(r)} 

= ’^kl,k 2 (?^), 


nr) 


( 20 ) 


where, 77*2 = (^ 1 ^ + kl^)/2 is the pair four-momenta of the quasiparticles in medium, 
and ■u^(r) = 7v[l,v(r)], /ii, 2 (r), and T(r) are the source four-velocity, the pair chemical 
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potential, and the source temperature at particle freeze-out, respectively. Equation (IT^ 
gives the relationship between k*^u^{r) and k^u^{r), which is used in calculating /ki,k 2 (?") 
and nki,k 2 (?^)- 


III. BBC RESULTS FOR HYDRODYNAMIC SOURCES 


A. Source distributions 


For hydrodynamic sources with a Bjorken cylinder, the four-dimension element of the 
freeze-out hypersurface can be written as 


^ //.('!', ri, J]) drifr^dr], 


( 21 ) 


where r, r^, and rj are the proper time, transverse coordinate, and space-time rapidity of the 
element. The function /^(r, r^, rj) is related to the freeze-out mechanism that is considered, 
and v) corresponds to the source distributions of proper time and space in the 

calculations [see Eqs. flT^ and flTT]) ]. In this work we assume that ^ mesons are frozen out 
at a fixed temperature Tf and use the AZHYDRO technique Q, H, Q to calculate the 
freeze-out hypersurface element. 
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FIG. 1: Distributions of the freeze-out points of 4> mesons in the z = 0 plane for the initial 
conditions eo = 8 and 20 GeV /ira? and Rx = Ry = 4: fm. 
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We show in Fig. [T]the distributions of the freeze-out points (source points) of 0 mesons in 
the z = 0 plane for the initial conditions cq = 8 and 20 GeV/fm^ and Rx = Ry = 4: fm. The 
distribution prohle for the lower Ty is wider than that for the higher Tf because of the source 
expansion. And the distribution prohles increase with increasing initial energy density. In 
Fig. [21 we show the normalized distributions of the transverse coordinate and time of the 0 
freeze-out points in the z = 0 plane, which are obtained by projecting the two-dimensional 
distributions in Fig. [1] to the coordinate and time axes, respectively. One can see that 
the transverse-coordinate distributions are similar to Gaussian distributions. The temporal 
distributions increase with time nonlinearly and have sharp falls to 0 at long evolving times. 
The widths of the spatial and temporal distributions increase with increasing initial energy 
density and decrease with increasing freeze-out temperature. 


Tf=165MeV Tf=140MeV 



Tf=165MeV Tf=140MeV 



FIG. 2: (Color online) Normalized distributions of transverse coordinate and time of the 0 freeze- 
out points in the z = 0 plane for the same initial conditions as in Fig. [H 


B. BBC functions 


In Fig. [3l we plot the BBG functions of 00 in k — m^ plane for the hydrodynamic sources 
with different initial radii and energy densities. In the calculations, we take p in the region 
(—1,1) and the chemical potential of boson-antiboson pairs, /ii, 2 (r) = 0. The freeze-out 
temperature of 0 meson is taken to be 140 MeV 0-7|. The variation of the BBG functions 
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(a) Rx =3 fm, Ry = 3 fm, eo=8 GeV/fm^ (b) R^ = 3 fm, Ry =4 fnti, eo=S GeV/frr? 



FIG. 3: BBC functions of (/>(/> for the hydrodynamic sources with different initial conditions and 
Tf = 140 MeV. 

with the modihed mass m* is similar to that of the BBC functions calculated previously 0- 
. However, the BBC functions for hydrodynamic sources exhibit oscillations as a function 
of k, compared to the BBC functions for sources with Gaussian spatial distributions and a 
temporal distribution of the exponential decay 

To examine the reason for the oscillations of BBC functions, we calculate the BBC 
functions for sources with spatial distributions obtained from the hydrodynamic freeze-out 
points, dN/d\r±\, [see Figs. |2][a) and[2][b)], and the parameterized temporal distribution. 


F(t) ^ (a + br^) 0{r - Tq) - r), 


( 22 ) 


where a, b, and Tmax are three parameters. As shown in Fig. 01 temporal distributions 


Fi and F 2 are similar to the distributions in Fig. [2][d). 


n Fig. 01 the thin solid line 
7||. It is much different from 


is for the temporal distribution of the exponential decay 
the distributions of Eq. 0221 ) . We plot in Fig. 0]the BBC functions of (jxj) for sources with 
hydrodynamic spatial distributions [see Fig. [2][b)] and parameterized temporal distributions 
of Eq. fl2^ and exponential decay (see Fig. 0]). Here, Figs. [5](a) and[5][b) are with temporal 
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6 L-fi,2 = («+^'t^^)e('t-'to)6('^max-T:) ^ 

; -fj(a=1.0, fo=0.035,Tujjj^=7) ; 

5r - f’2(«=0-9.^=0025,T„^=10) 2 

4 L -^■3=5exp[-(T-To)/2]e(T-To) J 



T (fm/c) 


FIG. 4: (Color online) Parameterized temporal distributions Fi and F 2 in Eq. (j22p for the param¬ 
eter sets [a = 1.0 (fm/c)“^, b = 0.035 (fm/c)“^, Tmax = 7 fm/c] and [a = 0.9 (fm/c)“^, b = 0.025 
(fm/c)“^, Tmax = 10 fm/cj. The thin solid (blue) line F 3 is for the temporal distribution of the 
exponential decay. 


(a) Rx =4 fm, Ry=4 fm, 60=8 GeV/fm^ 


(a') Rx=4 fm, Ry =4 fm, eo =8 GeV/fm^ 




(b) Rx=4 fm, Ry=4 fm, eo=20 GeV/fm^ 


(b') Rx=4 fm, Ry=4 fm, eo=20 GeV/fm^ 




FIG. 5: BBC functions of cjicj) for the sources with hydrodynamic spatial distributions [see Fig. 
[2][b)] and parameterized temporal distributions (see Fig. S]). (a, b) With temporal distributions 
Fi and F 2 in Fig. [H (a', b') with temporal distribution F 3 in Fig. 01 


distributions Fi and F 2 in Fig. 01 and Figs. [5](a') and[5](b') are with temporal distribution 
F 3 in Fig. 01 In the calculations, the source temperature is taken to be Tj = 140 MeV, 
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and the source velocities used are still the hydrodynamic source velocities at the freeze-out 
points. One can see that there are oscillations in the BBC functions calculated with temporal 
distributions Fi and F 2 , which have sharp falls to 0 at long times. However, the oscillations 
disappear in the BBC functions calculated with the temporal distribution of exponential 
decay, F^. The peak values of the BBC functions for temporal distributions Fi and F 2 are 
smaller than those of the BBC functions for temporal distribution T 3 . The reason, as will 
be seen, is that the width of temporal distribution T 3 is much smaller than those of Fi and 
F2. 


In the calculations of the BBC functions with parameterized temporal distributions, we 
have 


C(k, -k) = 1 + ^ 



J^DvWI‘(k)dr, 

2 


'lAs(- 

k) dr] 


where 


D,ik) = 




(23) 


(24) 


\V2-Vl\ lF(T)dT ’ 

^^(k) = j [4^_k(’^)c-k,k(r)nk,_k(r) + Ck,_k(r’)slk,k(^)K,-k(r) + l)]dh±, (25) 

I^(k) = y^^[|ck,k(r)pnk,k(’^) + |sk,k(^^)P(?^k,k(?^) + (26) 

Cki,k2(’^) = cosh[/ki,k2(^^)], Ski,k2(^^) = sinh[/ki,k2(^^)]- (27) 


For the parameterized temporal distribution in Eq. 022 p . we have 

[b^d^ -|- 4a(a -|- M^) sin^((i\/r? 7 ?”+~P coshr^)] 


1/2 


\D,(k)\ = 


(28) 


2\/m? -|- k‘^ cosh 77 {ad + bd^/3)\ri2 — rjil 
where d = ( xma v — uq), and the approximation b -C a{m? + k!^){cos\irjY/{chY is taken. The 
oscillations of the BBC functions are from Dr,{k). We plot \Dr^{k) \ in Fig. | 6 ]for the parame¬ 
terized temporal distributions in Eq. fl22p (thick lines) and the temporal distribution of expo¬ 
nential decay (thin lines), in this case \Dr^{k)\ = {(?72 — hi)^[l + 4(m^-|-fc^) cosher] Af^]} 
and At = 2 fm/c. We take rji = —1 and ri 2 = 1 in the calculations. The magnitudes of Dr^{k) 
behave as oscillations for the parameterized temporal distributions Fi because of the sharp 
falls in the distributions at long times (see thick lines in Fig. 0]). However, the magnitudes 
of D^{k) for the temporal distribution of the exponential decay, F 3 in Fig. 01 are smoothed. 
The magnitude of D^{k) for F 3 are larger than those for Fi because the width of the F^ 
distribution is smaller than that of Fi. For all the temporal distributions, the magnitude of 
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Djj{k) decreases with increasing rj. And, the BBC functions are determined by the product 
of Dr^(k) and I^{k). 



FIG. 6 : (Color online) Magnitudes of D^^k) for temporal distributions Fi and F 3 in Fig. 01 


C. Dependence of the BBC function on the direction of the particle momentum 

For anisotropic sources, the BBC functions depend not only on the magnitude of the 
particle momentum, but also on its direction [^. This is because the source velocities in 
different directions lead to the variation of with the momentum direction. To examine 
the dependence of the BBC functions on the direction of the particle momentum, we use 

cosa = fc^/|k|, cos/3 = fca;/|kr|, ^|kr| = ^k^ + k'^'^ , (29) 

to describe the direction of the particle momentum, a and /3 are the polar angle and 
azimuthal angle of the particle. In Fig. [3, we plot the BBC functions of </>(/) in the cos/3-/c 
plane for hydrodynamic sources with different initial radii and energy densities. Here, n\ 
is taken as 1.05 GeV, corresponding to approximately the peaks of the BBC functions (see 
Fig. [3]). For < Ry, the source expanding velocity in the x direction is higher than that 
in the y direction. The BBC functions for sources with R^ < Ry are larger at cos/3 = 0 
than those at cos/3 = 1, because the average values of are smaller at cos/3 = 0 

than that at cos /3 = 1 in this case [^. For the source with R^ = Ry, the BBC functions are 
independent of cos/3. 

In Fig. [HI we show the dependence of the average BBC functions, {C{k, —k))^, of (fxf) on 
the cosine of the particle azimuthal angle for hydrodynamic sources with different initial radii 
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(c) Rx =3 fm, Ry =4 fm, eo=20 GeV/fm^ 


(d) Rx =4 fm, Ry =4 fm, £o=20 GeV/fm^ 




FIG. 7: BBC functions of (pcf) in the cos/3-k plane for hydrodynamic sources with different initial 
radii and energy densities. Here, m* is taken as 1.05 GeV. 
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FIG. 8: (Color online) Dependence of the average BBC functions, {C{k, —k))k, of (fxj) on the cosine 
of the particle azimuthal angle for hydrodynamic sources with different initial radii and energy 
densities. Here, m* is taken as 1.05 GeV and the momentum region averaged is 0-1 GeV/c. 


and energy densities. Here, m* is taken as 1.05 GeV and the momentnm region averaged 
is 0-1 GeV/c. The BBG fnnctions are independent of the azimnthal angle for transverse 
isotropic sonrces. However, the BBG fnnctions for transverse anisotropic sonrces increase 
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with increasing azimnthal angle of the particles (0 < (3 < vr/2). For transverse isotropic 
sources, the average BBC function for the higher initial energy density is smaller than that 
for the lower initial energy density. The reasons are that the source with a higher initial 
energy density has a larger expansion velocity and a wider temporal distribution of freeze- 
out points. For transverse anisotropic sources, we observe that the average BBC function for 
the higher initial energy density is larger than that for the lower initial energy density. This 
is mainly because the many oscillations of the BBC function for the source with the higher 
initial energy density lead to an increase in the average value in the momentum region [see 
Figs. [71[a) and[71[c)]. 


(a) Rx=3 fm, Ry =4 fm, eo=8 GeV/fm^ (b) Rx =4 fm, Ry =4 fm, eo=8 GeV/fm^ 



FIG. 9: BBC functions of (/></> in the cos a-k plane for hydrodynamic sources with different initial 
radii and energy densities. Here, m* is taken as 1.05 GeV. 


We plot in Fig. [9] the BBC functions of (jycf) in the cos a-k plane for hydrodynamic sources 
with different initial radii and energy densities, and for n\ = 1.05 GeV. The ridge values 
of the BBC functions for a hxed k decrease with increasing cos a. The reasons for this are 
that the average source longitudinal velocity is higher than the average source transverse 
velocity for hydrodynamic sources with Bjorken longitudinal boost invariance and the 
higher longitudinal velocity leads to smaller average values of at cos a = 0 than at 
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cos a = 1 [7[. 

Because cos a is related to particle pseudorapidity by 

^ = tanh“^(cos a), (30) 

the polar angle dependence of the BBC functions can lead to pseudorapidity dependence 
of the BBC functions. In Fig. [101 we show the dependence of the average BBC functions, 
{C{k, —k))k, of 00 on the pseudorapidity of the particle for hydrodynamic sources with 
different initial radii and energy densities. Here, m* is taken as 1.05 GeV and the momentum 
region averaged is 0-1 GeV/c. The BBG functions decrease with increasing |^| as expected. 
Because the average transverse velocities for sources with a higher initial energy density 
are higher than those for sources with a lower initial energy density, the BBG functions for 
sources with a higher initial energy density decrease more slowly with increasing \y\. 
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FIG. 10: (Color online) Dependence of the average BBC functions {C{k, —k))k of (pcf) on the 
pseudorapidity of the particle for hydrodynamic sources with different initial radii and energy 
densities. Here, m* is taken as 1.05 GeV and the momentum region averaged is 0-1 GeV/c. 


IV. SUMMARY AND DISCUSSION 


In the hot and dense hadronic sources formed in high-energy heavy-ion collisions, particle 
interactions in medium might lead to a squeezed BBG of boson-antiboson pairs. The in¬ 


vestigations o 
distributions 


0the BBG in previous works are for sources with space-time-separated source 
7|. The smoothed temporal distribution of the exponential decay scales the 
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BBC functions and leads to monotonic BBC functions with respect to the particle momen¬ 
tum js-?]. 

Relativistic hydrodynamics is successful in describing the source space-time evolution in 
high-energy heavy-ion collisions. In this paper, we investigate the BBC functions of 00 
for sources evolving hydrodynamically in (2 -|- 1) dimensions and with longitudinal boost 
invariance. For hydrodynamic sources, the BBC functions oscillate as a function of the 
particle momentum. The reason for the oscillations is that the temporal distributions of 
sources evolving hydrodynamically have sharp falls to 0 at large evolving times, compared 
to the temporal distribution of the exponential decay. We also investigate the dependences 
of the BBC functions on the directions of the particle momentum. For transverse anisotropic 
sources, the anisotropic source velocity leads to the dependence of the BBC functions on 
the particle azimuthal angle. The BBC functions are minimum when the azimuthal angles 
of the particles reach 0. Because the average source longitudinal velocity is higher than 
the average source transverse velocity, the BBC functions increase with decreasing absolute 
value of the particle pseudorapidity. The oscillations and the dependences on the particle 
azimuthal angle and pseudorapidity are the signihcant signatures for detecting the BBC in 
high-energy heavy-ion collisions. 

For ki = k, k2 = — ki = — k, the BBC function (^(k, — k) for the 

hydrodynamic source is related to the temporal Fourier transformation of the space-time 
distribution of source freeze-out points [see Eqs. (jll), flT^ . and flMll ]. So, the BBC func¬ 
tion is very sensitive to the temporal distribution of the source, and an appropriate source 
space-time distribution is important for estimating the BBC effect in high-energy-heavy 
ion collisions. In Refs. {2 7|, a sudden freeze-out assumption at time Tf is adopted, and 
then a parameterized distribution of Tf (exponential decay) is used in the calculations to 
suppress the BBC functions. In Ref. [2^, the author argues the appropriateness of the 
sudden freeze-out assumption of time and discusses the exponential suppression of the BBC 
function based on three parameterized temporal distributions. In this work, we extract the 
space-time distributions of the source from the (2 -|- l)-dimensional hydrodynamics. The 
temporal and spatial distributions of the source are related in the calculations. Although 
it is a signihcant advance compared to the parameterized space-independent distributions 


of time used in 


M 


, further investigations based on more realistic models in which 


the model parameters are determined by the experimental data on the observables, such as 
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single particle spectra, elliptic flow, HBT radii, are needed for the expectations of the BBC 
effect in high-energy heavy-ion collisions. 
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